On low-temperature structural phase transitions 
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We comment on zero- and low-temperature structural phase transitions, expecting that these 
comments might be relevant not only for this structural case. We first consider a textbook model 
whose classical version is the only model for which the Landau theory of phase transitions and 
the concept of "soft mode" introduced by Ginzburg are exact. Within this model, we reveal the 
effects of quantum fluctuations and thermal ones at low temperatures. To do so, the knowledge of 
the dynamics of the model is needed. However, as already was emphasized by Ginzburg et al. in 
eighties, a realistic theory for such a dynamics at high temperatures is lacking, what also seems to be 
the case in the low temperature regime. Consequently, some theoretical conclusions turn out to be 
dependent on the assumptions on this dynamics. We illustrate this point with the low-temperature 
phase diagram, and discuss some unexpected shortcomings of the continuous medium approaches. 



I. INTRODUCTION 

Zero- and low-temperature (T) phase transitions are 
nowadays a subject of great interest (see, e.g., Refs. [1-5] 
for recent reviews). The special case of structural phase 
transitions deserves, in our opinion, a special attention. 
First, it is very convenient when introducing the topic 
of low-T phase transitions although, to the best of our 
knowledge, this pedagogical facet of the structural case 
has not been developed in the literature. One of the 
purposes of the present paper is just to develop this facet. 
Second, the discussion of structural transitions allows to 
reveal some unsolved problems which might have a fairly 
broad interest. 

It is worth mentioning that our study will be restricted 
to the region of small fluctuations (not very close to the 
phase-transition point). This region normally is not the 
region of main interest in the aforementioned papers, but 
the main specific features of the phase transition anoma- 
lies are clearly seen already in this region, not to mention 
that for interpretation of the experimental data this re- 
gion is quite often the most relevant one. 

A considerable part of the theory of low-T structural 
phase transitions is very simple. Its formulation uses ele- 
mentary formulas of quantum and statistical mechanics, 
and its development involves a fairly simple mathemat- 
ics. Nevertheless, this elementary theory suffices to dis- 
cuss some points of general interest such as the validity 
of the Landau theory, the soft-mode concept, the role 
of quantum fluctuations in defining the phase-transition 
point, the specific features of the low-T phase diagram, 
etc. This constitutes the first part of the paper where, be- 
cause of pedagogical considerations, we use a very simple 
model. Nevertheless, even within this elementary treat- 
ment, there arise some questions as well as not completely 
justified assumptions which will be discussed in the sec- 
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ond part of the paper. 

These questions and assumptions refer to the charac- 
ter of the dynamics of the order parameter near the zero- 
and low-T phase transitions. This character has not been 
successfully explained for high-T phase transitions: the 
origin of the so-called central peak in the soft mode spec- 
trum is understood only qualitatively [6]. For zero- and 
low-T structural transitions this question has not been 
studied at all, although the dynamics of the order pa- 
rameter is much more important here. Indeed, according 
to the classical statistical mechanics the static properties 
of the system do not depend on its dynamics. This is be- 
cause (gaussian) integration over momenta simply gives 
a factor in the corresponding partition function. But the 
situation is different when quantum effects play a role. 
In this case, the partition function does not factorize be- 
cause momenta and coordinates, now operators, do not 
commute with each other [7] . Therefore, a lack of exact 
knowledge of the dynamics impedes obtaining definite re- 
sults for, e.g., such a "static" property as dependence of 
the phase transition temperature on a control parame- 
ter (e.g., strain or pressure) in the low-T region. Given 
this situation, we will discuss several possibilities without 
proposing a finite conclusion about which of them corre- 
sponds to the reality. For this discussion we need no 
model at all, and the system is considered in this second 
part as a continuous medium. 



II. THE SINGLE-ION MODEL 

The so-called single- ion model (see, e.g., Ref. [8]) 
is very convenient when illustrating a zero-T structural 
phase transition. Within this model one assumes, first of 
all, that the crystal is composed by two types of atoms, 
say A and B. Our aim is to describe "active" A-atoms in 
the simplest way, so we further assume that i) the sub- 
lattice of i?-atoms can only be deformed homogeneously 
and ii) the interaction between A atoms is a nearest- 
neighbor interaction mediated by springs. Additionally, 
there is an interaction between A and B atoms which is 
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responsible for the relative position of the corresponding 
sublattices. Restricting ourselves to the orthorhombic 
case, let us choose the unit cell with J3-atoms placed at 
the apices of the corresponding cell (see Fig. 1). Thus, 
the potential acting on A-atoms due to the B ones has to 
be symmetric with respect to the center of this cell. This 
is so if this potential has i) a minimum in the center of 
the unit cell or ii) two symmetric out-of-center minima. 
In the following we shall assume that i) is the case when 
the crystal is strongly compressed and then, along the 
z-axis, it turns into case ii) with diminishing the com- 
pression (see Fig. 1). This makes possible a change in 
the mean position of A-atoms, i.e. a phase transition, in 
a fairly simple way. 

The potential energy of the system then can be written 

as 
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FIG. 1: The model: unit cell, effective potential acting on 
A-atoms due to B ones, and an illustration of the interaction 
between A-atoms. 
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where ur represents the displacement of the A-atom 
along the z-axis in the Jf?th unit cell. The first sum in this 
expression represents the effective potential acting on A- 
atoms due to B ones. Let us characterize the compression 
of the system by the magnitude w = (Vq — V)/Vq, where 
V is the volume of the system and Vq is this volume at 
zero pressure for the (noncquilibrium) configuration in 
which all A-atoms are maintained in the center of the 
corresponding unit cells (i.e., ur = 0). Thus, by taking 
a = a(w — wq), with a > 0, and b as a positive constant; 
Wq gives the strain at which the form of this potential 
change from one-well to two-well. [The usually small dif- 
ference between V and V (\w\, \w \ -C 1) turns out to be 
relevant for the change in the sign of the coefficient a only, 
so we shall not distinguish between V and Vq anywhere 
but here.] The second sum in Eq. (1) is the interaction 
potential between A-atoms, where c is the stiffness coeffi- 
cient of springs linking pairs of A-atoms (see Fig. 1) and 
summation is carried out over nearest-neighbors only. 



U = U + N[-u< + -ul 



(2) 



where N is the number unit cells. Minimizing this po- 
tential we find the equilibrium value of uq: 



2 JO (w > w ), 
"°' eq -\-f (w<w ), 

and the corresponding value of the potential energy: 
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In accordance with these formulas, the change in the form 
of the effective potential acting on A-atoms and the phase 
transition take place simultaneously at wq- 



Phase transition anomalies 



A. Static properties: A classical zero-T transition 

Let us suppose at this point that the mass of A-atoms 
is infinite, so they can be treated as classical particles. 
Consequently, the configuration of the system will be the 
one which simply minimizes the potential energy. The 
static properties of the system will be in accordance with 
this configuration, so let us proceed to determine it. 

It is clear that the minimum of the potential energy 
corresponds to the configuration in which the springs 
linking A-atoms do not experience any deformation. So 
all the atoms will be located in the same minimum of the 
effective potential created by B-atoms: ur = uq. Eq. 



As we have mentioned, the sublattice of £?-atoms, and 
therefore the system as a whole, can be compressed ho- 
mogeneously by applying pressure. Let us see how the 
corresponding stiffness of the system changes as a result 
of the phase transition. To this end, we present the po- 
tential energy of the system with A-atoms at the center 
of the corresponding unit cell as 

Uo ^U° + ^ + ... (5) 

(the term linear in w is absent by virtue of the defini- 
tion of w), and further take into account the energy of 
the mechanism applying the pressure: — P(Vq — V) = 
— PVqw. Minimizing the total energy with respect to w 



we find this pressure: 
Cow 



P = 



(w > w ), 



Cow - n^r(w - w ) (w<w ), 



(0) 



where Co = U /V and n = N/V. As a result, the stiffness 
is 



= dP_ = |Co ^ (w > wo), 
dw; |Co-nf^ (w<w ). 



(7) 



Its behavior is illustrated in Fig. 2, which is just the 
standard behavior for an anomaly described within the 
Landau theory of second-order phase transitions. 



2. Exactness of the Landau theory for the model 
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FIG. 2: 



having the meaning of the Fourier components of the 
displacement field: 



(9) 



It is worth noticing that we have no fluctuations within 
this model, neither thermal (as long as T = 0), nor quan- 
tum (as long as the masses are infinite). This makes 
Landau theory [8-10] to be an exact theory here. It is 
frequently said that the Landau theory is a mean-field 
approximation, but this example shows that this state- 
ment is highly inappropriate. 

The possibility of having a phase transition with no 
fluctuations questions another frequent saying: "zero- 
T phase transitions are due to quantum fluctuations." 
Why? We had no fluctuations but we had a transition. 
Equally unjustified is to say that high-T phase transi- 
tions are due to thermal fluctuations. Thermal fluctua- 
tions give rise to some contribution, of course, but as a 
rule they are not the unique reason of the transition. The 
origin these sayings may be in the fact that phase tran- 
sitions are frequently illustrated using the Ising model. 
For this model, it is true that the transition is due to 
the thermal fluctuations only. But while the Ising model 
played an important role in the theory of phase transi- 
tions, it is very specific, reflecting (as any model) no more 
than an aspect of much a more many-faceted reality. 



B. Dynamics 

Let us now consider that the masses of A-atoms are 
finite, although very large; so that the motion of these 
atoms is possible, but the corresponding dynamics can 
still be considered as classical. Let us study some specific 
features of such a dynamics. To do so, it is convenient to 
introduce the new variables 



In terms of these variables the potential energy Eq. (1) 
has the form 

k 

+ ^ U kU k >U k nU_ k _ W _ k n , (10) 

k,k',k" 

where a(k) = a + 4c(sin 2 %^ + sin 2 ^ + sin 2 %^), with 
I being the cell parameter (here it has been taken into 
account that u-k = u k ). By separating the zero Fourier 
component of the displacement field: uq = N^ 1 J2r u R> 
which evidently has the meaning of the mean value of the 
displacement field, the above expression for the potential 
energy can be rewritten as 

U = U Q + N(lul + \ui) + NY j a -^\ukf 

fc#0 

+ Nbu 2J UfcUfc'W-fc-fc' 

k,k'=i0 

+ ^ UfcUfc'Wfc"U-fc-fc'-fc"- (11) 

k.k' \k"^0 

where a(k, u ) = a(k) + 36uq. Below we shall consider 
the case of small-amplitude oscillations such that the con- 
tribution of the last two sums in this expression can be 
neglected. It is worth noticing that this does not mean 
to neglect all the anharmonicity of the system as long as 
the coefficient b is still present in the remaining terms: 

U^Uo + N^ul + ^+N^-^lukf. (12) 

fc#0 
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R 



(8) 



In the simplest case, the equation of motion for the 
Fourier components of the displacement field can be writ- 
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of harmonic oscillators [see Eq. (13)]. For a fixed value 
of Mo, the ground state energy of the system then can be 
written as 



E = U + 



Twj c {k, uq) 



(16) 



(recall that T = 0). As a result of the possible displace- 
ments of A-atoms along the z-axis, here we have the 
classical (macroscopic) contribution considered in pre- 
vious section plus a quantum term accounting for the 
ground-state energy of the corresponding oscillators with 
normal frequencies ui c (k,uo). In further calculations it 
is convenient to take the continuous-medium limit of 
model (i.e., to replace summation by integration 
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FIG. 3: 
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ten as 



The phase-transition point 



mu k + a(k, u )uk = 0, 



(13) 



where m is the mass of the A-atoms. We then have 
an optical branch with the dispersion law uj 2 (k,Uo) = 
a(k,uo)/m (see Fig. 3). For small wavevectors: 



w^(fc, Uq) — (a + 3bu 2 + ck 2 )/m, 



(14) 



where c = cl 2 . Substituting here the equilibrium value 
of uo, Eq. (3), we obtain the normal frequencies of the 
system as a function of the control parameter w. For 
k = that is 



Wc(0,U O ,eq) 



' a\w — w \ / \ 
2ol\w— wq\ / , \ 

—4s — ■ ( w < w o)- 



(15) 



This behavior is illustrated in Fig. 3. It is worth noticing 
that, at w = wq, this optic branch has the same disper- 
sion law as the acoustic one: oj oc k. As a result of this 
behavior, first noticed by Ginzburg [11], this k = mode 
is usually termed as the soft mode associated with the 
transition. 



C. Zero-T transition: Quantum effects 

Let us now consider that the masses of A-atoms are 
"normal" (not specially large). In this case, if we consider 
A-atoms as immobile atoms when studying the thermo- 
dynamic properties of the system we are wrong. Even at 
T = these atoms are moving exhibiting what is called 
quantum fluctuations, and we have to take into account 
this quantum effects. 

In a first step, the system can be considered as a set 



Because of the dependence on uq of the quantum term 
(oc K) in Eq. (17), it can be said that there effectively 
is a quantum renormalization of the function E(uq). By 
expanding the last term in Eq. (17) in power series of 
Mo we can obtain, for example, the quantum correction 
to the coefficient a. The resulting coefficient, a*, is of 
particular interest to us: when a* = the system losses 
its stability with respect to nonvanishing values of Mo, 
what just defines the phase-transition point. Taking into 
account that 



duj c (k,u ) 
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u =0 



d 2 w c {k,uo) 
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u =0 



muj c (k, 0) 



we find that 



dk 



3hbv 



2m J (2tt) 3 w c (fc,0)' 



(18a) 
(18b) 

(19) 



where v = V/N = I 3 is the volume of the unit cell. When 
trying to determine the phase-transition point one real- 
izes that this formula is, however, somewhat contradic- 
tory. The coefficient a has to be negative in order to 
further obtain a* = 0. But in this case, in accordance 
with Eq. (14), there is a region of wavevectors for which 

LU 2 < 0. 

This difficulty can be circumvented if this region is 
much smaller than the whole Brillouin zone. When this 
possibility takes place the system can be labeled as dis- 
placive and, when calculating the integral in Eq. (19), 
one can put a = quite safely. As a result one finds that 
the phase-transition strain is given by w c — wo + 5w, 
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where Sw is the contribution due to quantum effects 
(zero-point quantum fluctuations) which can be esti- 
mated as 



Sw 



3hbv 



k 2 dk 



3hb 



(20) 



Alternatively, by realizing that the frequency of the soft 
mode should vanish at the "experimental" value of phase 
transition strain w c , one could replace a — > a* in the inte- 
gral in Eq. (19). This is equivalent to take into account, 
in an effective way, higher order corrections in Eq. (19) 
(see below). 

Although the problem seems to be overcame, it is in- 
structive to discuss this "overcoming" in more detail. 
What we have done is to "correct" the dynamical moduli 
of the system in a meaningful way. This is quite similar 
to what is done when calculating, e.g., the phase tran- 
sition temperature of a classical displacivc system: here 
they are the static moduli what are corrected in a proper 
way. By realizing that our system is in fact an anhar- 
monic system, as indicated by the last two sums in Eq. 
(11), it can be expected that going beyond the approxi- 
mation of decoupled oscillators such a correction (renor- 
malization) will be obtained. The possibility of using 
a perturbation theory to account for this anharmonicity 
is therefore quite attractive, and valuable irrespective of 
the character, classical or quantum, of the problem in 
question. However, the quantum case is somewhat more 
complex than the classic one as long as dynamics plays 
a key role. One can realize that we did two things in 
fact: i) assumed that dynamics is still well-described in 
the same terms (i.e., within a soft-mode scenario) and 
then ii) corrected the corresponding dynamical moduli. 
The former has not to be necessarily true in a real case. 
Indeed, one of the well known shortcomings of the the- 
ory of high-T structural phase transitions is that the soft 
mode behavior has never been observed experimentally 
to a full extent: the "soft mode", at best, diminishes its 
frequency considerably; but this frequency never goes to 
zero at the phase transition point. At the same time, 
close to the phase transition point there appears the so- 
called "central peak" in the order parameter fluctuation 
spectrum, and the increase of the order parameter fluctu- 
ations turns out to be comprised within this peak. Such 
a peak is a natural feature of the so-called order-disorder 
systems, for which the order-parameter dynamics is re- 
laxational instead of phonon-like as we have considered 
so far. Though there is no consistent analytical theory 
of the central peak, its appearance is not surprising at 
all [6] . Any system with a phase transition is inherently 
anharmonic. So the phonon-like ("soft mode") picture is 
not the full picture; there is another part, essentially an 
anharmonic one, which reveals itself in the central peak. 
Similar effects are quite possible at low Ts: the presen- 
tation of the motion as a set of normal vibrations, even 
followed by taking into account the anharmonism within 
a perturbation theory, may fail to describe accurately 



the dynamics of the system close to the corresponding 
phase-transition point. Consequently, one can do noth- 
ing but assume Eq. (19) with caution and hope that it 
is qualitatively correct. 

The importance of the quantum corrections to phase- 
transition strain can be estimated as follows. The very 
possibility of observing the transition implies that \wo \ is 
much less than the atomic value of this magnitude |u> a t | ~ 
1. The other constants in the model, however, can have 
their "atomic" values. Introducing the atomic (binding) 
energy e at ~ h 2 /(m e l 2 ), where m e is the electron mass 
and I ~ lA the atomic distance, we then have 



a 



£at 



Eat 



Thus wc find that 



Sw 


^ /m e \V 2 




Wq 


V m J 


w 



(21) 



(22) 



As long as m e <C m, but |u> a t| 3> |u>o|, here we see that 
the correction to the classical phase-transition strain can 
be quite important: it may give rise to a magnitude w c 
comparable, or even greater, than Wq. 



2. Order-disorder (spin-like) limit 



We have discussed the displacive limit but what hap- 
pens if the value of \a\ at the phase-transition point is so 
large that lo 2 < in most of the Brillouin zone? In this 
case, one simply has to realize that the starting point 
was not chosen appropriately. Looking back we can see 
that this case corresponds to an effective potential acting 
on A-atoms due to B ones with two profound wells even 
in the symmetric phase. Under these circumstances, in- 
stead of characterizing the motion of ^4-atoms by their 
displacement from the center of the corresponding unit 
cell (what ur really means), it makes more sense to as- 
sume, first of all, the atoms to be confined inside one of 
the two wells and then take into account the possibility of 
delocallization by, e.g., quantum tunneling. This can be 
done by associating the two initial positions of ^4-atoms 
with the two possible orientations of (pseudo-) spins 1/2 
and accounting for the presence of a transversal field in 
order to reproduce the tunneling. Thus our model be- 
comes a spin model exhibiting a zero-T phase transition 
(Ref. [2, 4]). The corresponding excitations (spin waves) 
are similar to the optical vibrations considered above in 
the sense that there is a "soft spin- wave" whose frequency 
goes to zero at the phase-transition point. In this sense, 
when it comes to zero-T transitions, displacive and order- 
disorder limiting cases are not so different one another. 
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3. Phase-transition anomalies 

Let us now calculate the quantum contribution to the 
anomaly in the stiffness of the system within the dis- 
placive scenario. The main purpose of such a calculation 
is to estimate the region in which the Landau theory is 
applicable, so it suffices to consider the symmetric phase 
(uo = 0). From Eqs. (16) and (14), we find that 



c = 



1 d 2 F 
Co 



Tina 2 



ha 2 



dk 



1 



In 



327r 2 m 1 / 2 c 3 / 2 \a(w — wq) 



47T 2 C 



3/2 



(fc,0) 



(23) 



where £o is in accordance with Eq. (7). Let us analyze 
this result in some detail. 

At first sight, it seems that the stiffness of the sys- 
tem changes its sign and diverges at w = wq. However, 
as long as we are considering the lowest order correc- 
tion only, we simply can say that it diminishes close to 
wq. In any case, why close to wq? An anomaly at the 
phase-transition point would be not surprising but, in 
accordance with Eq. (19), the phase-transition point is 
not Wq- What happens? The answer is that here have 
the same kind of inconsistency that we have found before 
[compare the integral in Eqs. (19) and (23)]. So the re- 
placement of wo by the actual value associated with the 
phase-transition strain, which can be extracted from the 
experiments, solves it. 

An estimate of the region of applicability of the Landau 
theory in our case can be made by demanding that the 
(quantum) contribution to the stiffness of the system is 
much smaller than the discontinuity obtained within the 
Landau theory itself. That is 



Tina 



32 7 r 2 m 1 / 2 c 3 / 2 



In 



Air 2 c 



a(w — w c ) 



< n 



2b 



(24) 



[see Eq. (23)]. In accordance with previous estimates 
of the coefficients appearing in this expression [see Eq. 
(21)], this gives 



In 



W at 



< ltr ( — 

m e 



1/2 



10 4 



(25) 



As wee see, it is hardly possible to abandon this region 
in real experiments. 

It is worth mentioning that the behavior of the stiffness 
of our model in the scaling region has been known since 
long ago [12]. Roughly speaking, it is such that one has 
to replace In by In 1 / 3 in Eq. (23). As wc see, both the 
first-order approximation and the theory for the scaling 
region provide qualitatively the same results (see Fig. 4). 
However, as a rule, the first-order approximation gives 
rise to an overestimate of the relevant magnitudes in the 
very vicinity of the phase-transition point: in our case, by 
extrapolating the increase of the stiffness obtained within 



the first-order approximation to the scaling region one 
overestimates such an increase. 



4- Thermal effects 

Thermal effects can be revealed by following the same 
considerations as above. At finite temperature, the free 
energy of the system can be written as 



fa o b , 
F = F + N[-u 2 + -u 4 Q 



V 



dk hLo c (k,u ) 
(2tt) 3 2 



TV 



dk 



ln{l-exp[-ftw c (fc,ii )/T]}. (26) 



(Notice that this expression is nothing but the free en- 
ergy of a set of harmonic oscillators, see, e.g., Ref. [9].) 
Therefore, the coefficient at the term quadratic in uq is 



STibv 



dk 



2m J (2tt) 3 w c (fe,0) 
2>Tibv f dk n[u> c (k, 0)] 
1W (2tt) 3 w c (fe,0) ; 



(27) 



where n(ui) = [exp(huj/T) — 1] _1 is the Bose-Einstein 
distribution function. Here we have the quantum contri- 
bution computed in previous section [see Eq. (19)] plus 
a thermal one (the latter integral). Similar to what we 
have done when discussing the quantum contribution to 
W c , it is reasonable to replace a — > a* in the expression 
of ijj c entering in Eq. (19). 

Let us discuss the phase diagram. The border between 
the symmetric and non-symmetric phases in the (T, w\ 
plane is defined by the condition a* = 0. This gives the 
line 



w - w c = -a V (T c ) , 



(28) 



in the (T, w)-plane, where / (T c ) is the last integral in 
Eq. (27) evaluated at a* = 0. At low temperatures 
[Ulu(0,Q) «T« fuj(k max , 0)], this integral can be esti- 
mated as 



dk n[w c (fe,0)] If n[w c (k,0)] 2 



(2tt) 3 w c (fc,0) 2tt 2 7 u) c (k,0) 



1 /.kmaxCT) 



2tt 2 J hcu 2 (k,0) 



k 2 dk 



,3/2 



2-k 2 Ti 2 <?I 2 v 



(29) 



taking into account that the Bose-Einstein distribution 
function decreases very strongly within the integration 

1/2 

interval [it vanishes for k > fc m ax(T 1 ) = Prza T\. This 
gives 



w c oc T 



(30) 



At high temperatures [hu> c (k, 0) <C T], the Bose- 
Einstein distribution can safely be replaced by its classi- 
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imcnts the central-peak phenomenon appears time and 
again, and this lacks a theoretical explanation beyond 
"hand- waving" arguments [6]. The situation for zero- 
and low-T phase transitions is unknown, but one may 
suspect that it is similar. In consequence it is quite rea- 
sonable to consider, within the continuous medium ap- 
proach, both hypothetical cases: that of the soft-mode 
(phonon-like dynamics) scenario and the central-peak- 
like (relaxational dynamics) one. As we shall see, the 
aforementioned universality is lost in the second case for 
which the continuous medium approach becomes no more 
than a model. 



cal limit: 

dk n[oj c (k,0)] 



(2tt)3 W c (M) 

We then have 



k 2 dk T 



2ir 2 TuJ 2 (k,0) nhev 



W — W r OC T r 



T. 

(31) 

(32) 



III. CONTINUOUS MEDIUM APPROACH 

At this point, it is worth making a comparison be- 
tween origin of the quantum and thermal contributions 
in a* obtained previously [see Eqs. (19) and (27)]. While 
the quantum contribution comes from the ground-state 
energy of all the optical phonons, the thermal contri- 
bution is obtained from the optical phonons with small 
wavevectors only. In consequence, the actual value of 
phase-transition strain at zero temperature (w c ) is sen- 
sitive to the microscopic details of the system, whereas 
the function T c (w) defining the phase transition temper- 
ature as a function of strain is not. This can be taken 
as a justification to use a continuous medium approach 
when studying thermal dependencies. Just in this sense 
one speaks about an universal low-T behavior of the sys- 
tems, exemplified by the temperature dependence of spe- 
cific heat of solids (the Debye law): it comes from the 
contribution of small-wavevector acoustic phonons and, 
therefore, is characterized by macroscopic quantities such 
as the velocity of sound. 

However, in our case, this universality is not so ev- 
ident. It can be expected when the soft-mode scenario 
takes place, as we were able to argue for a classical zero-T 
transition. But as soon as we make the model more re- 
alistic, taking into account real values of the ion masses 
and, therefore, the quantum effects, as well as consid- 
ering other degrees of freedom, in particular, acoustic 
modes, the dynamics of the system becomes far more 
complicated. This has not been studied to a full extent, 
for which there is a reasonable explanation: difficulties 
similar to those found for the high-T order-parameter dy- 
namics arc indeed expected. As we have mentioned, the 
soft-mode scenario is not completely applicable to high- 
T phase transitions: in both real and numerical cxper- 



A. Phonon-like dynamics 

Our arguments about the universality were referred 
just for the phonon-like case. In this case, the the- 
ory can be made somewhat more realistic, even model- 
independent, than the one presented above where only 
one ("active") optical branch is taken into account. It 
is straightforward, in particular, to take into account the 
acoustic branches (e.g., to account for the motion of B- 
atoms as well). Indeed one may expect that these acous- 
tic phonons give rise to a significant contribution to a*: 
they are low-frequency excitations of the system. Con- 
sider the free energy 



,'a*(T) , b , 



dk huj ac (k,u ) 
(2tt) 3 2 



3TV 



dk 



ln{l - exp[-7kj ac (fc,u )/T]} 



(33) 



where a*(T) is given by Eq. (27). This is nothing but Eq. 
(26) where the contribution of the acoustic phonons has 
been added explicitly. This contribution depends on uq 
because of the corresponding dependence of the velocity 
of sound: 



w ac (/c,u ) = c(u )k, 



(34) 



which can be taken as c(u ) = c (l + ^Ug). (This velocity 
has to be understood as an averaged velocity of sound in 
the same sense as in Ref. [9].) Thus the new coefficient 
of the free energy at the quadratic term in uq can be 
written as 



a**{T) 



a*(T) 
3hav 



2tt 2 



- + n[cj ac (fc, 0)] | w ac (fc, 0)k 2 dk. 

(35) 



As we see in this expression, we have a new contribution 
to the phase-transition strain w c which is not surprising 
(the part of the integral which does not vanish at T = 0). 
But wc also have a thermal one which can be estimated 
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as 

(36) 

Having in mind that the thermal activation of the optical 
phonons obtained above is exponentially suppressed with 
moving away the phase-transition point [see Eq. (27)], 
one realizes that this thermal contribution due to acous- 
tic phonons may be the most important one in a sig- 
nificant region of the (T, w)-plane. However, this is not 
the case close enough to the phase transition point where 
a *(T) — a*(0) oc T 2 . Consequently, the form of the phase 
diagram does not change from the one expounded above. 

These changes take place if, for instance, there ex- 
ist long-range interactions in the system. In this case, 
complicated dispersion relations may arise such that the 
dimensionality of the integrals in Eq. (27) effectively 
increase. In the case of uniaxial ferroelectrics, for in- 
stance, dipolar interactions lead to uj 2 (k, 0) = (a + ck 2 + 
4ir cos 2 9)/m for small wavevectors. This further gives 
a*(T) - a*(0) cx T 3 (sec, e.g., Rcfs. [3, 13]). The lin- 
ear coupling between the order parameter and one of the 
component of the strain tensor, i.e., the piezoelectric ef- 
fect in the case of ferroelectrics, gives rise a peculiar de- 
pendence a*(T) - a*(0) cx T 5 / 2 [14]. These long-range 
interactions, however, do not modify the high-T behav- 
ior. 



B. Dissipative dynamics 

Let us now consider that there is a central peak in 
the order-parameter fluctuation spectrum, at least close 
to the phase-transition point. Such an assumption is far 
from being trivial: it implies that there is an essential fre- 
quency dispersion in the order-parameter response func- 
tion and, therefore, a dissipative dynamics of the order 



parameter even at T = 0. In a number of the systems, 
e.g., for perfect crystals without phase transitions, this 
is certainly not the case. But the situation with such 
an inherently anharmonic system as the one we are con- 
sidering here (recall that a < in our case) is unclear. 
What seems to be clear, however, is that for crystals with 
some defects the dissipative dynamics of the type we are 
supposing here takes place even at T = [15]. Therefore, 
even if our consideration would prove to be irrelevant for 
pure systems, it might still be relevant for systems with 
some defects. A dissipative dynamics for a set of vari- 
ables means that these variables are not coordinates of 
a closed system but there is a "reservoir" where the sys- 
tem in question transfers its energy to. In our case, the 
long-wave acoustic phonons could play the role of such a 
reservoir. When calculating the contribution of the set of 
dissipative variables one should, of course, not forget the 
contribution of the reservoir. Only in the case when the 
latter contribution is irrelevant it makes sense to consider 
a system of dissipative variables as the main contributor 
(see, e.g., Ref. [16]). Let us mention that even apart from 
possibility of the central peak, an account of a phonon 
damping at T — due to e.g. defects is of some interest 
for zero- and low-T transitions. We shall begin with this 
case by postulating the equation of motion 

mu k + 7u fc + a(fc, u )u k = 0, (37) 

where the viscosity coefficient 7 accounts for the dissipa- 
tion. 



1. Thermal effects 

Making use of some results obtained by previous au- 
thors (see, e.g., Ref. [16]), we can study the case in which 
the motion of the "active" -atoms is governed by Eq. (37). 
In this case, the free energy can be written as [16] 



F = F 



TV 



dk 

J2nf 



In 



/27T[X 1 (k,uo)X 2 (k,u )} 1 / 2 



lnT 1 



Ai(fc,it ) 



lnT 1 



A 2 (fc,M ) N 

v J 
(38) 



where v = 2nT/Ti and {Xi(k,uo)} are the roots of the 
equation A 2 (fc, u ) — jX(k, u ) +uj 2 (k, u n ) = 0, satisfying 
the relations 

Ai(fc,u ) + A 2 (fe,u ) = 7, (39a) 
Ai(fc, u )X 2 {k 7 u Q ) = uj 2 (k, u ). (39b) 

It is worth mentioning that the first term Fq in Eq. (38) 
now includes the contribution due to the reservoir. In 



principle, this contribution could exhibit a nontrivial be- 
havior close to the transition point as a result of the 
coupling between the active atoms (which clearly "feel" 
the transition) and the degrees of freedom in the reser- 
voir (which also "feel" the transition because of this cou- 
pling). Nevertheless, if these degrees of freedom forming 
the reservoir are low-frequency acoustic phonons, which 
would be quite natural, it can be shown that their contri- 
bution i) is not substantially modified because of damp- 
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ing [15] and ii) is not the most important one close to 
the (low-T) phase-transition point as we have seen in 
Sec. Ill A [14]. 

The above formulas allow us to compute the relevant 
magnitudes close to the phase transition by following the 
same procedure as before (see Ref. [14]). The coeffi- 
cient at the quadratic term in uq in the free energy, for 
instance, is found to be 



a*{T) 



Zbv 



dk 



+ 



mk B T J (2tt) 3 
^[l + Ai(fc,0)/i/] 



1 



2 Wc 2(fc,0) 

^[l + A a (fc,0)/i/] 



i/[Ai(fe,0)- A a (M)] 



(40) 



Close to the phase-transition point, this further gives 

a*(T) — a*(0) oc T 2 , (41) 

in the relaxation limit (m — > 0). It is worth mention- 
ing that this behavior is found irrespective of long-range 
forces [14], which is related to the fact that the thermal 
activation of the elementary excitations of the system is 
substantially modified as a result of the damping [no- 
tice that the Bose-Einstein distribution function in Eq. 
(27) has been replaced by the psi functions in Eq. (40)]. 
In fact, the most important conclusion here is that, as 
a result of damping, the paradigm: low-T behaviors = 
small wavevectors, is not longer valid. This is because all 
degrees of freedom may give a contribution to the corre- 



sponding thermal behavior as a result of the damping [if 
m — > 0, in Eq. (40) there is no temperature-dependent 
cutoff analogous to that in Eq. (27)]. In some sense, 
this is similar to what happens in the case of the phase 
transition strain: one has to go beyond the continuous 
media approach to calculate this strain. The situation is 
new, however, in the sense that this happens even for the 
thermal behavior. 



IV. CONCLUSIONS 

The subject of low-temperature phase transitions has 
been examined, dealing with the structural case and fo- 
cusing on the region of small fluctuations. Within an 
oversimplified but illustrative model we have discussed 
some points of general interest, such as the validity of 
the Landau theory, the role of quantum fluctuations in 
defining the phase-transition point and the specific fea- 
tures of the low-temperature phase diagram. We have 
pointed out also that a profound study of dynamics, both 
experimental and theoretical, is still needed to get fully 
justified results about the phase diagrams and anomalies 
at structural and may be other low-temperature phase 
transitions. 
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